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TECHNICAL PAPER 


SOLUTION OF THE ANGLES-ONLY SATELLITE TRACKING PROBLEM 

INTRODUCTION 


The determination of the orbit of a satellite vehicle or ballistic missile is obviously of great 
practical importance. In the usual case, the azimuth, elevation, range, and range-rate of the target are 
measured from an Earth-fixed radar site. When range and range -rate are not available, the determination 
of the orbit from the azimuth(s), elevation(s), and the time(s) of the observation is referred to as the 
“angles-only problem.” 

A good summary of this problem is treated in reference 1 . In that source, it is pointed out that the 
problem was solved by Gauss, Laplace, and Escobal. These solutions seem to be far better known 
among astronomers than among aerospace engineers.* 

In the present report, another solution is derived which is particularly well suited to numerical 
calculations. It is shown that all of the orbital elements (semi-major axis (a), eccentricity (e), argument 
of perigee (S), inclination (7), and longitude of the ascending node (SI)) can be derived from five 
measurements of the azimuth (and the times of observation), or five measurements of the elevation (and 
the times of observation). Once these elements have been isolated, the epoch can be derived from a 
given time measurement. 

An advantage of the present formulation is that it can be readily extended to the case of an 
observer aboard a second satellite rather than an Earth-fixed observer. Since only line-of-sight angles 
would be required to determine the orbit of a second satellite, the technique could be used to eliminate 
range and range-rate hardware that is presently required for rendezvous. This option will be treated in a 
future paper. 


THE DRIVER PROGRAM 


In the case of an actual radar installation, the azimuth and elevation angles, as well as the time of 
the observations, would be recorded once a target is acquired. Since this report deals with theoretical 
sightings rather than actual observations, it is necessary to build a driver program to simulate what a real 
radar station would observe. To that end, a vehicle having arbitrary orbital elements is assumed and the 
details of the sightings from a general configuration are derived. From a general configuration, specific 
orbital elements are used for simulation. Once the elevation, azimuth, and time are calculated, the orbital 
parameters are hidden from the isolation scheme. The isolation then has to rederive a complete 
specification of the orbit from the table of elevation angle versus time (or the azimuth versus time). 
Since a ballistic missile that impacts the Earth is obviously in a “terminated” orbit, the same theory 
applies in either case, i.e., one can track an incoming missile using the present theory. For the ballistic 
missile case, drag would be a significant factor, however. Drag is not included in the present report. 


* The author wishes to thank Dr. John Hanson, EL-58, MSFC, for pointing out this literature on this subject. 



Let (?, y, k) describe an Earth-centered inertial coordinate system and locate a tracking station by 
the polar coordinates ( R a , <f) 0 , 6 0 ) (fig. 1). We explicitly recognize the Earth’s rotation by setting 
(j) 0 = At the tracking station, a spherical coordinate system is constructed via the unit vectors 
(fi r , uq, uq). Note that this coordinate system is modified from the usual convention in that the u e points 
due south along the meridian rather than the conventional practice of orienting it to the north. The target 
is located, relative to the tracking station by a distance, p, and two angles, the azimuth a, and the 
elevation angle, /3. 

The unit vector transformation is: 

u r - cos 0 o ■ cos 1 + cos 6 0 -sin p 0 y+sin 0 o ■ k , 

“0 = sin (p 0 •?+ cos (j) 0 -y , 

Uq = sin 0 O -cos (p Q -i+sin 6 0 -sin (j> 0 -y'-cos 6 0 ■ k , 
along with the vector from the center of the Earth to the radar station, R 0 , as 

R 0 = R 0 cos 0 o cos <f) 0 1+R C ;Cos 6 0 sin (f) 0 y+^sin 6 0 • k , (1) 


where R 0 = | R () |. 

In the local coordinates of the radar station, the p vector is given by 

p = p (sin (j -ft r +cos sin a iL+cos ft cos <X-uq), (2) 


where p = | p |. 

Inserting the expressions for u r , , and u e into the expressions for p , one finds 

p = p - [(sin P -cos 6 0 -cos (f> 0 -cos P -sin a -sin (p a -i-cos P -cos a -sin 6 0 -cos (j) 0 )■? 

+ (sin P -cos 0 O -sin (j) 0 -i-cos P -sin a -cos q ) 0 +cos p -cos a -sin 6 0 -sin (j> 0 )-y 

+(sin P sin 6 0 -cos p cos a cos 0 O )-k] . (3) 

Define a vector r from the center of the Earth to the satellite. For reference 2, one can write r in 
terms of the inclination, /, the right ascension of the ascending node, SI, the argument of perigee, w, and 
the true anomaly,/, as 

r = x-i+y ■ j+z-k = r-[cos (/+w)-cos <$-sin (/+© )-sin <^-cos I] i 
+r-[cos (/+©)-sin ^2+sin (/+ <w)-cos <^-cos 7] y 

-i-r-sin (/+Si)-sin I k , (4) 

where r = \ r \ (see fig. 2). Also, one can write 

r = a-(l-e 2 )/(l+e-cosf) = p/(l+e-cosf ) , (5) 
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where a is the semi-major axis of the orbit and e as the eccentricity and p = a- (1-e 2 ). 

Vectorially, 

? = R 0 +p. (6) 

Equating the components of equation (6) from equations (5), (3), and (1), 

p -(sin p cos 0 O 'cos (p 0 -cos /3 -sin a -sin (f> 0 -i-cos ft -cos a -sin 6 0 -cos (f> 0 ) 

= r-[cos (/+ S))-cos c^-sin (/+®)-sin $-cos 7]-7? 0 -cos 6 0 -cos (j) 0 , (7) 

p (sin p cos 6 0 sin +cos p sin a cos (j) Q +COS p cos a -sin 6 0 -sin <j ) 0 ) 

= r-[cos (/+S)-sin $+sin (/+2))-cos <$-cos 7]-7? 0 - cos 6 0 -sin <j> 0 , (8) 

p (sin /?-sin 0 o -cos /fcos a-cos 6 0 ) = r-sin (/+2})-sin 7-7? 0 -sin 6 0 . (9) 

Taking the square root of the sum of the squares of equations (7), (8), and (9) gives 
p = (r 2 +/J^-2-r-7? 0 -{[cos (/+S)-cos 0 o -cos (<$-$, )-sin (/+S)] [cos 0 O -cos 7 sin (&l-<p 0 ) 

-sin 7-sin 6 0 ]}) 1/2 . (10) 

Although the preceding equations can be solved directly for a and p, a more subtle approach 
yields additional quadrant information that can be valuable during numerical calculations. To this end, 
one first calculates (8) -cos 0 o -(7)-sin (j) 0 to get 

p-cos /?-sin a= r-[cos (/+ 5))-sin (<$-$, )+sin (/+ 5i) cos 7-cos (<$-^ )] . (11) 

Next, calculate (7)-cos <p a +( 8)-sin (p 0 to yield 

p (sin P -cos d 0 +cos P -cos a -sin 0 o ) 

=r-[cos (/+S)-cos )-sin (/+©)• cos 7-sin (< Q-<j> 0 )-R a cos 6 0 ] . (12) 

From (12) -sin 0 o -(9)-cos 6 0 one has 

p-cos p - cos a- r-{cos (/+S))-cos (c )-sin G c 
-sin (/+£&)• [cos 7-sin (<$-$, )-sin 6 0 +sin 7-cos 0 o ]} . (13) 

Now, divide equation (1 1) by equation (13) and obtain the formula for the tangent of a as 

tan a= [cos (/+2)) sin (<$-$, )+sin (/+S)-cos 7-cos {Sl~<p 0 )]/{cos (/+S)-cos (Q-fy )-sin 6 0 

-sin (/+2))-[cos 7-sin (<Q-^ )-sin 6 0 +sin 7-cos 6 0 ]} . (14) 

Equation (14) may be rewritten in many forms, with one of the more useful being 
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tan a - [sin (<$-0, )+tan (/+ 2))-cos /-cos (Sl-<p 0 )]/{cos 6 0 - {cos (c $l~<p 0 )-tan 6 0 

-tan (/+8))-[cos /-sin )-tan 6 0 +sin /]}) . (15) 

In equation (14), one can separately interrogate the numerator and denominator (interpreted as 
sin a and cos a ) to obtain the needed quadrant information on a 4 

The elevation, (5, can be obtained in several ways. One such method would come from 
multiplying equation (12) by cos 6 0 and equation (9) by sin 6 0 and adding to obtain 

p-sin / = r-{cos (/+5>)- cos ($-$, )-cos 6 0 -sin (/+8))-[cos 7-sin (6 2-0, )-cos 6 0 

-sin 7-sin 6 0 ]}~Ro • (16) 

Squaring and adding equations (11) and (13) yields a formula for cos (5 , but this process destroys 
information. Dividing equation (15) by that equation for cos [3 will yield a tangent formula which 
involves r and R 0 but not p. 

Another approach to determine ft is simply to utilize equation (10) to determine p and then use 
either equation (11) or equation (13) to obtain / as an arc cosine. Equation (11) would be greatly 
preferred since cos a will usually vanish during a satellite pass. 

The magnitude of the radius vector, r, which appears in many of the above equations involves 
the semi-major axis, a, and the eccentricity, e. If one has full knowledge of the orbital elements a, e, 7, 
to, and SI, and if one also knows the epoch (or a related time point), one can now obtain the values of p, 
a, and (5. The time reference serves to determine the true anomaly,/, and the position of the observation 
site, since 0 O = 0 o (O- 

For given values of the orbital elements and epoch (needed for r and 0 O ), the path of computation 
is now clear. Equation (10) yields the values of p, and equation (15) can be solved for a. Once p and a 
are known, equation (11) provides the value of /. 

It is apparent that not all values of the epoch will yield values of (p, a, (3) which correspond to 
physical data that could be recorded at a given observation site. Subsidiary conditions must be imposed 
on the epoch (as well as on 0 O ) in order to guarantee that the satellite is above the horizon of a ground 
plane which is tangent to the site. (“Over the horizon radar” is not considered here, but such a device 
could easily be incorporated into these equations.) Physically, apparition of the target is expected when 
[3 = 0 and occultation when [3 = n. These conditions will be fulfilled if 

p-R o =0. 


Dotting equation (1) with equation (3) gives 

cos (/+5>)-cos (Sl-(j > 0 )- cos 0 o -sin (/+5))-[cos 7-sin ($-p 0 )-cos Q 0 - sin 7 -sin 0 O ] =RJr. (17) 

This equation is time dependent through the variables/, r, and <j> 0 . In order to recognize the explicit time 
dependence of the longitude of the observation site, write 


0 O 0 O +0-^- 


t The author wishes to thank Dr. Larry Mullins, MSFC, NASA, for pointing out the existence of this formula. 
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( 18 ) 



Since the true anomaly, /, is time dependent, r must be expanded by equation (5). In terms of the 
eccentric anomaly, £, 


r= a- [1-e-cos (£)]. ( 19 ) 

The eccentric anomaly is related to the time by Gauss’ equation (zero referenced to the time of perigee 
passage) as 

t = • ( £— <? sin £) f (20) 

and to the true anomaly,/, via 

/ = 2-tan- 1 {/[(l+e)/(l-e)]-[tan (£ 12 )] } . (21) 

Since we are in the driver program, it is assumed that the orbital elements are all known. One 
must isolate a value of the time such that equation (17) is satisfied, along with equations (18) through 
(21). Since no useful analytic solution exists even for equation (20), there is little hope of solving 
equation (17) analytically. For this reason, a Newton iteration was developed that defined 

H = cos (/+6>) cos (Sl~(f) o )-cos 6 0 -sin (/+©)• [cos /-sin (c$-^ )-cos 6 0 

-sin / sin d 0 \-RJr, (22) 

in accord with equation (17). To isolate a zero of H, the iterator requires that one also know dH/dt. 

dH/dt = -{sin(/+6))cos ( Q-fy ) cos 0 o +cos (/+©)• [cos / sin (6l~<p 0 ) cos 6 a 

- sin /-sin 0 o ]}■/(/?•//) /r 2 -{ cos (/+S)-sin (<$-$, )-cos 6 0 

+sin (/+£y)-[cos /-cos (<^-0 o )-cos d 0 ]}-^+{RJp)-e - sin f-J(p-n)/r 2 , (23) 

where 

df/dt = \J\p-a-(\-e 2 )]/r 2 = *J(p-p)/r 2 , 
and 

d/dt(\/r) = -(e/p )- sin f->J(p-/lj/r 2 

have been used. 

Although the use of equations (22) and (23) can isolate the apparition and occultation of the 
target, it is convenient to employ approximate calculations to estimate those times as starting values for 
the Newton iteration. This can be conveniently done if one assumes (here) that <f> is zero, but allow /to 
vary in accord with Keplerian dynamics. 

To approximate the apparition and occultation of the target, from equation (13), one defines 

X 0 = cos (Sl-(p 0 >cos e 0 , (24) 

X l = [cos I sin (Sl~<f> 0 )-cos 0 O -sin I - sin 0 O ] , (25) 

so that equation (13) can be written as 


X 0 -cos (/+ty)-X]-sin (f+co)-R 0 /r = 0. 


(26) 
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Expanding the double angle formulas in equation (22), and abbreviating 


(27) 

(28) 


X 2 = X 0 -cos TO-Xj-sin w , 

X 3 = X 0 -sin tn+X^cos w , 

allows equation (22) to be written as 

X 2 -cos /-X 3 -sin f-R 0 lr = 0 . (29) 

Inserting r from equation (5) into equation (29) and cross multiplying by the denominator of r 

gives 


cos fiX 2 -R 0 -e/p)~ sin f-X 3 = R a e/p. 

(30) 

Defining the variable T) by 



R 

= tan _1 [X 3 /(X 2 -R 0 -e/pJ], 

(31) 

and £ as 



£ =cos _1 

(R 0 -e/{p^[X 3 2 +(X 2 -R 0 -e/p) 2 ] }>, 

(32) 

allows one to write equation (30) as 

COS (/+ T] ) = COS £ . 

(33) 

Equation (29) has two solutions, either 

/=£“* h 

(34) 

or 

/ = Ijc-ig+rj). 

(35) 


Either the solution of equation (34) or the solution of equation (35) will correspond to apparition 
on a nonrotating Earth; the other will correspond to occultation. Although Gauss’ equation assumes that 
the initial time point occurs at perigee, it is not difficult to introduce a false time reference of zero at 
apparition by defining this point as t D and then writing all time values, t, as t-t 0 . Once the time of 
apparition has been isolated, this time can be used as a starting point for the equations since the 
calculation starts at that instant. The second solution for/ (which could be the solution of either equation 
(34) or (35)) then supplies a starting value for the Newton iteration of equation (22) so that the method 
will converge on occultation rather than apparition. 

An orbital table can now be produced by choosing an arbitrary number of points between the/, 
(apparition) and the/- (occultation). Since each /and the orbital elements are known, one can specify all 
the orbital parameters at these points. Table 1 shows the results of such a procedure for 15 points on a 
typical orbit. The orbital elements were arbitrarily taken as a = 7,420 km, e = 0.1, w = 25°, I = 35°, and 
= 33°; (j) 0 was taken as 77° and 9 0 as 28.5°, and the Earth’s rotation rate, <p , was 2^/(24-3,600). 
Apparition occurred at f 0 = 8.481° (122.153 seconds past perigee), and occultation occurred at 
ff- 50.855° (749.288 seconds past perigee). The remaining entries of table 1 were calculated by the use 
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of equations (5) for r, (10) for p, (15) for a, and (11) for £ can most easily be obtained from equation 
(19). The equations that were used to compute the values of the orbital angles (j) and 6 are from standard 
celestial mechanics 


and 


(j) = tan 1 { [tan (/+&>) -cos 7+tan $]/[l-tan (/+®) cos /-tan <$] } , 
6= tan -1 [tan /-sin (</>-<$)] • 


(36) 

(37) 


It is interesting to determine how large f} can become in the field of view. This can be done by 
requiring that 3 /?/3/ = 0, and then solving for the corresponding value of /. When this value of / is 
inserted into equation (11), the maximum value of /3 will be found. If one uses equation (11), that 
procedure is algebraically rather messy. It is better to proceed by abbreviating 

X = cos (/+S)) cos (Sl-<b )-cos 6 0 -sin (/+©)• [cos /-sin (<$-<& ) -cos 0 O 

-sin I sin 0 o ] , 

so that equations (10) and (11) can be written, respectively as 

p = sJ(r 2 +R 2 0 -2-r-R 0 -X) 

sin p = ( r-X-R 0 )/p = {r-X-R 0 )!sJ (r 2 +R 2 0 -l-rR 0 -X) 


(38) 


(39) 


which immediately yields 

tan p = (r-X-R 0 )/[r-\/ (1-X 2 ) ] 


(41) 


If we now take the derivative of (41), equate 3)3/3/ to zero and clear, we find that the requirement for /3 
to be a maximum is 


R 0 (l -X 2 )(dr/df)+r- [r-R 0 -X]-(dX/df) = 0. 


(42) 


From equation (15), we have 

dr/df = r 2 e - sin flp, (43) 

while equation (34) gives 

dX/df = -sin (/+ 2))-cos )-cos 0 O -cos (/+©)• [cos /-sin (<$-^ )-cos d a 

-sin /-sin 0 o ] . (44) 

Taking (34) cos (/+5))-(44) sin (/+5>) gives an equation which can be solved for dX/df as 

dX/df = X-cot (/+S)-cos (&l~<p 0 )-cos 6 0 -esc (f+a>) . (45) 

Inserting equations (15), (43), and (45) into (42) gives, finally. 


(46) 


X 2 -[e-cos <y+cos (/+ty)]-X-[cos (<$-$ )-cos 6 0 -(1+e-cos f)+(p/R 0 )-cos ( f+co )] 
+[(p/P 0 )-cos ( 6l-(j) 0 )-cos 6 0 +e-sin /-sin (/+©)] = 0. 
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Equation (46) is cubic in the trigonometric functions of /. Even if one substitutes 
sin 2 /— » 1-cos 2 /, the residual sin/terms will require an additional squaring operation, so the best that 
we can hope for is a sixth order equation in cos/. This virtually precludes an analytical solution for / but 
numerical iteration works very well. For a numerical iteration, it is easier to use equation (42) (along 
with equations (38), (43), and (44)) rather than equation (46) (along with equation (38)). Once/ has been 
isolated, the maximum value of p can be obtained from equation (41). 

A far easier question is to determine whether or not P = nil in the viewing field. From equation 
(41), p can be nil only if X = +1. Using equation (38), this requirement implies that 

cos (/+2>)-cos 6 0 -cos ($-0 o )-sin (/+£&)•[ cos 6 0 -cos 7-sin (£l~<p 0 ) 


- sin 7-sin 6 0 \ = ± 1 . (47) 

Setting 

tan A= [cos 6 0 -cos 7-sin ( &l~<p 0 )-sin 7-sin 6 0 ]/[cos 6 0 -cos (c$-^ )] 

= cos 7-tan )-sin 7-tan 6 a -sec ( Sl~(j) 0 ) , (48) 

then, using spherical trigonometry, 

tan A= cos 7-tan ($-$ )+sin 7-tan 7-tan (c$-^ ) = tan (c$-^ )-sec 7 . (49) 

Equation (47) can now be written as 

cos (/+S))-cos A-sin (/+5))-sin A = ±1 , 
or 

cos (f+m+X) = ±1 . (50) 

Thus, p will achieve an angle of n/2 (or 3 k 12) if either of the true anomalies 

/ = -tan -1 [tan (#-$ )-sec I]-cb , (51) 

or 

/ = 7T-tan -1 [tan (c )-sec 7]— S> , (52) 

fall between apparition and occultation. 

From equation (14), an additional datum is available. It is interesting to determine the value of 
the true anomaly, / which results in a = 0. For this condition, equation (14) requires that 

sin (#-0 o )+tan (/+5))-cos 7-cos (c$-^ ) = 0 , 
or 

tan (/+©) = -sec 7-tan (&l~(f) 0 ) . (53) 

If p is now taken to be the longitude of the satellite (as opposed (j) 0 , the longitude of the tracking station) 
a standard equation from celestial mechanics gives 

tan (/+©) = -sec 7-tan (Sl-p) . (54) 
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For equations (52) and (53) to hold simultaneously, one must then have <j> = <j) 0 as the condition for 
a = 0. At this time 


/ = -{tan 1 [sec /-tan )]+<!)} . (55) 

If the satellite was observed exactly at a = 0 (and the associated time), it would be possible to use this 
datum to eliminate one orbital element in terms of the remaining four. However, nonlinearities probably 
make such a procedure more trouble than it is worth. 


THE SOLUTION OF THE ANGLES-ONLY PROBLEM 


In order to obtain the data in table 1, it was assumed that the orbital elements were known, but 
that assumption was used only to simulate data that could be recorded by an observer with full 
knowledge of the orbit. A naive observer who is located at a radar station and restricted to measuring 
only elevation and azimuth angles (as well as the time of the observations) could record only the sort of 
information shown in table 2. The “angles only” problem can now be stated quantitatively as follows: 
given the five sets of readings in table 2 , isolate the five orbital elements, a, e, a, I, and SI that could 
have produced those readings. It will be shown that all of the data of that table are not needed. Indeed, 
one can isolate the required orbital elements from the five reading of either a and the time or j 8 and the 
time. But listing a and (5 and the time is redundant. In practice, since all of the reading will be noisy, it 
would be prudent to solve the problem twice, with one solution involving a and the other solution 
involving (5 and then to compare the two answers. Although only five readings can be used to produce 
the results, these five readings are mathematically precise and thus are “safe.” In actual field work, one 
would probably take many more readings and compute the elements of the orbits from many sets of five, 
then statistically validate the final answers. These considerations will not be pursued in the present 
report. 


To begin, assume any convenient set of values for a, e, do, /, and SI . It is required only that they 
be reasonable. In most cases, some guess can be made about the semi-major axis (for a satellite, require 
thata (l+e) > R 0 ), the eccentricity (say 0 < e < 1) and the inclination (if the observation station is at 
9 0 ,I>e o ). 

One could proceed in the calculations in the following way. Once a set of orbital elements have 
been assumed, one can obtain a value of £ from equation (20) for a given time. Equation (21) then 
produces the corresponding value off. Once /has been calculated, obtain p from equation (10) and a 
from equation (15). Equation (11) then gives the value of /3. If the guesses about the orbital elements are 
exact, then the recorded tabular values of a and /3 will agree with the calculated values of these variables 
and the problem is solved. 

Another procedure was found to be more convenient, however. Given a starting set of guesses 
(a, e, co, I, SI), first calculate Xq, ..., X 3 from equations (26) through (28). Additionally, calculate the 
following convenient quantities: 


X 4 = sin {Sl-<p 0 ) , 

(56) 

X 5 = cos / cos ( Sl~(p 0 ) , 

(57) 

X 6 = X 4 -cos S)+X 5 -sin 0 ) , 

(58) 


9 


X 7 = X 4 sin S)-X 5 cos S) , (59) 

X 8 = X 6 -tan /? -X 2 -sin a , (60) 

X 9 = X 7 -tan /? -X 3 -sin a , (61) 

X 10 = X s -(R 0 /p) sin ae , (62) 

X\i = \J (X 9 2 +X l0 2 ) , (63) 

where the tabulated values of a are utilized. In the general case [sin (ft) 0], obtain £ from: 

t, - tan -1 (X 10 /X 9 ) , (64) 


and rj from 


Tj - sin 1 [(/? 0 /p)sin a/X n ] . (65) 

The value of/ which corresponds to the given ft must be either 

f = Z-Tl, (66) 

or 

/ = 2k-(Z+71). (67) 

(The decision of which /to use requires some programming logic, governed by continuity.) Once the 
value of/is obtained, from the inverse of equation (21), 

£ = 2-tan" 1 { v /[(l-e)/(l+e)]-tan(//2)] } , (68) 

with equation (20) then supply the value of time which corresponds to the chosen set of orbital elements 
at the isolated value of true anomaly. 

Unless the orbital elements were guessed precisely, the time which has been computed will not 
agree with the time which was recorded as a corresponding value to each given a, (5, or the (a,/)) pair. 
What is needed is an iteration scheme which will modify the values of the orbital elements in a 
systematic manner until the calculated times agree with the actual times recorded at the observation site 
(or the times given by the driver program). Notationally, let T T$ 0 be the “true” times and let 
T$ be the times which are calculated from equation (16). Denote (/■ as 

Ui = Ti(a,e,(d,Sl,l)-Ti 0 (i= 0, ...,4), (69) 

and let an “i” subscript denote corresponding time-varying quantities (/, £, etc.). 

Suppose that the orbital elements (a, e, a>, I, SI) differ from the true values of the orbital 
elements by amounts A a, Ae, AS, A I, and A SI (respectively). If these deviations were known precisely, 
one would have 


U i(a+Aci, e+Ae, gh-Aq), Sl+ASl, I+Al) — 0. 


(70) 
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The obvious next step is to expand the £/,• set via Taylor’s theory to obtain a Newton iterator: 

Ui(a+Aa, e+Ae, I+AJ, Sh-AS, $+A$) = Uj(a, e,I,a>,Sl)+ (d U i/da) -Aa+(dUi/de) 

■Ae+(dU i /dI)-AI+(dU i /dm)Adj+(dU i /dSl)ASl+ ... = 0, (71) 

where (/ = 0, ..., 4). In matrix notation, with yj as any member of the ordered set of orbital elements 
deviations (A a, Ae, A I, AS, A SI), and xj and member of the ordered set of orbital elements {a, e, /, S, 
<Q), equation (71) becomes: 

(Uf = m/dxjHyf (i = 0, ..., 4), (j = 0, ..., 4) . (72) 

Unless | ( dUj/dxj ) | = 0, one can solve for 

(y j ) T = (dU i /dx j )- l -(Uf. (73) 

The values of a, e, /, w, and SI can now be modified by the substitutions 

a — > a+Aa 
e — » e+Ae 

I — > I+Al > (74) 

W -» W+Aw 

Sl^Sl+ASl > 

If the new values of the orbital elements do not drive all of the [/,• values to zero, the process can 
be continued until convergence is achieved. 

One still needs to obtain the derivatives of each of the (/, with respect to each of the orbital 
elements, i.e., the matrix elements of dUi/dxj. Temporarily suppressing the subscript, it is convenient 
to proceed through a “chain calculation” as: 

dU/da = (3/2)-VT«W-(£-e-sin £)+J(a 3 /tf (l-e •cos £)d£da 


= V («///) • [(3/2) (£-e svs\ £+a { 1-e-cos £) d£da], ( 75 ) 

dU/de = s/ (a 3 /ju) -[(1-ecos £)-d£de], (76) 

dU/da b = \J (a 3 //X) -[(1-ecos £) d£da>], ( 77 ) 

dU/dQ. = (a 3 / /l) -[(1 -e-cos £)-d£dSl] , (78) 

dU/dl = \/ (a 3 ///) - [(l-e- cos £)-d£dI] . ( 79 ) 
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From equation (68): 


d£/da = ( r/a)-df/da 

(80) 

d£/de = [r/ \J (1-e 2 ) ][(df/de)/a- sin f/p] ? 

(81) 

d£/da) = ( r/a)df/dcb , 

(82) 

d&dQ = ( r/a)-dfldSl , 

(83) 

3#37 = (r/a)-a//a/. 

(84) 


One will also need the derivatives of £ from equation (64) and T] from equation (65): 

dr] _ -[7? 0 -sin a] (a dX n /da+X n ) 

^ " [a-X u V(fl 2 -(l-e 2 ) 2 -Xfi-[/? 0 -sin (ct)] 2 )] ’ 


dr] _ { [R 0 - sin cc\/[X x v j\-e 2 )\Y{2-e-X n -(\-e 2 )dX n lde] 




v /(a 2 -(l- e 2 ) 2 -X2 ] _[/ ?( .s in («)] 2 )]) 


(86) 


3t7/0x = ([7/sin a ]/Xj O-dXj x fox)lsl (a 2 -( \-e 2 f-X n 2 -[R 0 ■sin a] 2 )], 
where x = 7, w, or <Q. For the derivatives of £, one has 

dfydx = [X 9 - 0X , 0 /3x)-X , 0 - (3X 9 /3x)]/[X 9 2 +X , 0 2 ] 

= [X 9 - (3X , 0 /3x)-X ] 0 - (3X 9 /3 jc)]/X ! , 2 ’ 

(for x any orbital element). 

Using the given expression of /, equation (66), one is now in a position to calculate the 
derivatives of/ from, say: 

df/dx = 3£ /dx-drj/dx . (89) 

The next task is to calculate the derivatives for Xq, ..., Xjj. These follow (with * as a generic 
orbital element): 


BX 0 /dSl = -sin ($-<& )’Cos 0 o , 

(90) 

dXg/dx = 0 , (x ^ c$), 

(91) 

aXj/3/ = -sin 7-sin (<$-<& )-cos ^ -cos 7-sin 6 0 , 

(92) 

dX x !dSl = cos 7-cos(<$-$, )-cos ^ , 

(93) 


(87) 


( 88 ) 
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9Xj/3x = 0, (x = a, e, or 23), (94) 

dX 2 /dx = (3X 0 /3x)-cos a>-(dX j/3x)-sin 23, (x ? 23), (95) 

3X 2 /32 3 = (dX 0 /3tt)cos 23-(dXj/32))sin 23-X 0 sin 23-Xj-cos 2 3 , (96) 

dX 3 /dx = (3X (} /3x)-sin 2>-(3X 1 /9x)-cos 2), (x ^ 23) , (97) 

dX 3 /dcb = (5Xo/223) sin 23+(3X 1 /323>cos 23 +X 0 -cos 23-Xj-sin 23, (98) 

3X 4 /3$ = cos ($-$,), (99) 

3X 4 /3x = 0, (x*$), (100) 

dX 5 /dI = -sin /-cos ($-$, ) , (101) 

3X 5 /0c(2 = -cos /-sin ($-$, ) , (102) 

dX 5 /dx = 0 (x ^ /, x ^ $), (103) 

dX 6 /dx = (3X 4 /3x)-cos 23+(3X 5 /3x)-sin 23, (x ^ 23) , (104) 

3X 6 /323 = (3X 4 /323) cos 23+(3X 5 /323) sin 2>-X 4 sin 23+X 5 cos 23 , (105) 

dX-j/dx = (3X 4 /3x)-sin 23-(3X 5 /3x)-cos 23, (x & 23), (106) 

dXj/dd) = (3X 4 /323)sin 23-(3X 5 /323)cos 23+X 4 cos 23+X 5 sin 23 , (107) 

dX 8 /dx = (dX 6 /dx)-tan (3 -(3X 2 /3x)-sin a , (108) 

dX 9 /dx = (3X 7 /3x)-tan /J-(3X 3 /3x)-sin a, (109) 

3X, (/3x = dX 8 /dx, (x^e), (110) 

dX {( Jde = 3X 8 /de+R 0 - sin a, (111) 

ax n /3x = [x 9 -(3x 9 /ax)+x 10 -(ax 10 /3x)]/x n . (i 12 ) 


If the program utilizes both the a and (3 table, the development is now complete. Once the orbital 
elements have been chosen, we can obtain the derivatives of Xo,..., Xu with respect to (a, e, 23 , 1, $) 
each orbital element from equations (90) through (112). Given these derivatives, use equations (85) 
through (88) to obtain the derivatives of 7) and £ with respect to the orbital elements, and subsequently, 
the derivatives of/ with respect to the same variables. That allows a computation of the derivatives of e 
(equations (80) though (84)) and finally the derivatives of U from equations (75) through (79). This 
yields the derivatives of a five-dimensional vector, say ( dUo/dx ). The process is repeated for each of the 
other data sets (a„ 7)), to construct the matrix shown in equation (72). Matrix manipulation then 

yields the corrections to the orbital elements from equation(s) (74). 
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If the information on the azimuth, a (or (3), is not used, the situation is more complicated since a 
must be taken as another variable and equations (60), (61) and (62) must account for the fact that a is 
defined in terms of the assumed orbital elements. From equation (11), one can use equations (58) and 


(59) to write 

sin a = r-(Z 6 - cos f-X-j-sin f)/(p -cos /? ), (113) 

while equation (10) becomes 

P = \J [r 2 +R„-2-r-R 0 -(X 2 -cos f-X 3 - sin /)] . (114) 

Equations (107) through (111) would be modified to read: 

aXg/djc = (dX 6 /dx) tan [) -(dX 2 /dx)-sm a-X 2 - cos a -da/dx , (115) 

dXg/dx = (9X 7 /9x)-tan (3 -(dX 3 /dx)-sin a-X 3 - cos a da/dx, (116) 

dX] { fdx = dX s /dx+R a l p )- cos a -da/dx (x^e), (117) 

9X, (Jde = dX^/de +R 0 sin a+(R 0 /p ) -cos a -da/dx (x = e). (118) 

From equation (114) comes 

dp/dx = {(r-R 0 )-dr/dx-r-R 0 -[ cos f-(dX 2 /dx-Xydf/d x) 

- sin f-(dX 3 /dx+X 2 -df/dx\ }/p, (119) 

which is needed to evaluate daJdx\ from equation (113) one obtains both 


cos a -da/dx = {dr/dx-(X 6 -cos f-X-j-sm f)+r-[(dX(/dx-X- 7 df/dx)-cos f 

-(X 6 -df/dx+dX 7 /dx)- sin /] }/(p -cos f5 )-r-(X 6 - cos /'-X 7 -sin f)-dp /dtl(p 2 'Cos (3 ) , (120) 

and 

cos a =\J [p 2 cos 2 ) 3 -r 2 -(X 6 -c os /-X 7 -sin /) 2 ] Kp ■ cos [ 3 ) . (121) 

In order to take the derivatives of a with respect to each orbital element, one would also need the 
derivatives of r with respect to the same variables. Using equation (19), one finds 

dr/da = 1-c-cos £+a-e-sin £-d£/da , (122) 

9 ride = a-cos £+a-e-sin £-d£/de . (123) 

Using equations (115) through (123) allows the measured values for a to be ignored, but 
ignoring the measured values is paid for with the need for increased computation. Similar comments are 
applicable to the case where tabular values of a are assumed to be known but /3 is treated as a 
computable quantity. In the numerical results that are presented in the next section, tabular values of 
both a and [3 were assumed to be known. 
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NUMERICAL RESULTS 


The above theory was programmed in COMMON LISP on a Symbolics 6370. That platform 
allows for rapid prototyping and ease of equation development. The entire set of equations performed 
quite well; the Newton iteration scheme can converge a correct set of orbital elements from rather bad 
guesses. 

There has been an implicit assumption that there is a unique set of orbital elements which will 
produce the given table of observations. Although no attempt was made to establish uniqueness, one 
case was inadvertently found where the solution was not mathematically unique, but the extraneous 
solution did not correspond to a physically realistic orbit. Specifically, for the inputs of a , (3, and time 
from a two-dimensional suborbital case (with true a = 6,018.739, e = 0.724, 2) = 225, 1 = 0, and Si = 0), 
started from a rather precise initial set of guesses that assumed a = 6,020, e = 0.5, a> = 350, 1=0, SI = 0. 
The convergence criteria was met with a = 9,666.27930, e = -0.17980, a> = -25.54499 , 1 = 0, and SI = 0. 
This is not a physically realistic orbit (e < 0, a > 0), however, and is easily eliminated by requiring that 
e > 0. Additional constraints should be added to the iteration scheme, i.e., require that a-{ 1+e) > R 0 , 
-nil <l< nil, and 0 < S[ < n. In the case of tracking a comet or an asteroid (which are usually elliptical 
with respect to the Sun but hyperbolic with respect to the Earth), alternative constraints obviously would 
have to be used. 

In regard to two-dimensional iterations, in general, note that the convergence scheme is valid for 
those cases provided that the reduced dimension of the iteration space is accounted for. A three- 
dimensional iteration, faced with a two-dimensional problem, will encounter a singular matrix unless 
precautions are taken. Furthermore, when a two-dimensional problem is undertaken, one must take 
a = 3n/2 due to the choice of the coordinate system. 

During the computation procedure, the equation for i), equation (31), was particularly 
treacherous in that it is rather easy for (R 0 /p)- sin a/X 11 to underflow or overflow from the range 
-1 < sin -1 < 1 and produce complex values for a. The LISP language easily handled the complex values, 
but it balked at attempting to double-float the results. This was handled by the crude, though effective, 
artifice of simply setting a to either 3n/2 or n/2 depending upon the direction in which the insult 
occurred. As the iteration scheme approached more accurate values of the orbital elements, the problem 
disappeared. Iteration to isolate a satellite orbit tends to monotone convergence, but isolation of very 
short range missiles (-200 km) becomes extremely difficult since the trajectories have eccentricities that 
approach unity. For those cases, the Newton iteration, as presented here, is virtually inapplicable. 

Another difficulty, often present in Newtonian iteration schemes, was that the corrections for A a, 
Ae, AS) , A I and A SI were often so large that the iteration scheme began to wander aimlessly. Two 
methods were used to circumvent this problem. One option measured the magnitude of the errors and 
adjusted a “creep factor” according to the size of the errors. Another scheme placed absolute bounds on 
the allowable changes that could be made by any variable in a given iteration. This report was written to 
demonstrate a new approach to the angles-only problem, so no attempt was made to develop an optimal 
iteration scheme. Since the angles-only problem is important in the real world, it is expected that specific 
iteration schemes to optimize convergence may be forthcoming in the future. 

Two illustrative examples are given here to demonstrate the use of the equations given earlier. As 
mentioned previously, the first of these is a rather standard satellite orbit which is numerically described 
in Table 1; the fact that 16 sets of data are recorded is purely arbitrary. A subset of the data displayed in 
Table 1 is shown in Table 2. This is the entire set of data sent to the iteration scheme. As explained, the 
numerical examples used the values of a and (3 along the time of the observations, so that neither a nor 
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ft were calculated from the other as a part of the study. The iteration scheme converged the proper 
elements from the input guesses in 26 iterations (Table 3). Note that the figure of merit for the 
convergence (last column, Table 3) is monotone decreasing. The error calculation is defined as 

V U$+uf+uZ+l%+U} . 

Table 4 demonstrates the same data as Table 1, but this time for a suborbital vehicle. The case 
was deliberately taken as equatorial so that a two-dimensional iteration could be demonstrated. Since 
only three orbital elements (a, e, cb ) need to be isolated, one can make-do with only three readings. For 
a two-dimensional case, if the tracking station is located in the plane of the orbit, the azimuth angle, a, is 
meaningless. Thus, only three readings are needed. Again, the end points were /? = 0 or ft = n were 
included, but this is arbitrary. Table 4 shows the results of the iteration scheme from very poor initial 
guesses; indeed, the guesses are about as bad as they can meaningfully be. The iteration scheme does not 
display monotone convergence and appears to wander aimless for 61 of the 107 iterations that were 
needed to converge the case. 

An attempt was made to further stress the iteration scheme by requesting an isolation of the 200- 
km range suborbital planar case, but convergence was never attained. Certainly, the minimal Newton 
iteration presented here can be strengthened by any of several known methods. 


CONCLUSIONS 


The preceding theory has derived another solution to the “angles-only” problem. This method, 
possibly used in conjunction with the alternative methods covered in reference 1, has the potential to 
decrease the amount of hardware that is currently required for rendezvous and missile tracking. 

The numerical work that is presented demonstrates that the theory yields a practical scheme that 
actually solves the problem. Unfortunately, it also indicates that additional research is required to extend 
numerical convergence to eccentricities which are very close to unity. 

A planned future paper will relax the requirement of an Earth-fixed observer and treat the case of 
an observer in an orbiting satellite. 
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Table 1. First example. 

SEMI-fIRJOR AXIS (A)=7420. KM, ECC=. 10000 INC( I ) =35 . 000 DEG., RRG . PERIGEE(u)=25.000 DEG., LONG. RSC . NODE (0) =33 . 000 DEG. 
APPARITION OCCURS AT f=8.4Bl deg. [T=122.154 sec.] RND OCCULTATION OCCURS AT f=50.309 deg. [T=748.883 sec.] 


TRUE RNOMOLY, DEG. 

RADIUS, KM. 

RANGE, KM. 

ECC. RNOMOLY, DEG. 

TIME, SEC. 

0 f =0 . 48133600 

r=6684. 64603564 

P=2001. 40179057 

6=7.67420293 

T=122. 15382475 

1 fell. 33937028 

r=6689 . B7187027 

P=1702. 55541086 

6=10.26292799 

T=163. 40242366 

2 f =14.19230929 

r=6696. 58128091 

p=1406. 84597518 

6=12.84934444 

T=204. 65102256 

3 f =17 . 0309138B 

r=6704. 75468903 

P=1117 . 35905908 

6=15.43288791 

T=245. 89962147 

4 f =19. 87797554 

r=6714. 36839048 

P=841. 02426446 

6=18.01300611 

T=287. 14822037 

5 f =22. 708321 98 

r=6725. 39470844 

p=596. 64903616 

6=20.58916104 

1=328.39681928 

6 f =25.528B2210 

r=6737. 80216907 

P=440. 75876006 

6=23.16083097 

1=369.64541818 

7 f =28. 33839058 

r=6751. 55569755 

P=470. 56098302 

6=25.72751237 

T=410. 89401708 

8 f =31. 13599196 

r=6766. 61683174 

P=660. 49278407 

6=28.28872158 

T=452. 14261599 

9 f=33. 92064421 

r=6782. 94395090 

P=914. 90251289 

6=30.84399635 

T=493. 39121489 

10 f=36. 69142178 

r=6800. 49251644 

p=1192. 33738686 

6=33.39289717 

T=534. 63981380 

11 f=39. 44745811 

r=6819. 21532193 

P=1479. 03113160 

6=35.93500846 

1=575.88841270 

12 f=42. 18794762 

r=6839. 06274951 

p=1769. 64478692 

6=38.46993949 

1=617.13701161 

13 f=44. 91214716 

r=6859. 98302993 

P=2061. 70019529 

6=40.99732514 

T=658. 38561051 

14 f=47. 61937693 

r=6801. 92250357 

P=2353. 86307735 

6=43.51682652 

T=699. 63420942 

15 f =50. 30902095 

r=6904. 82588003 

p=2645. 32323647 

6=46.02813134 

T=740. 88280832 

AZIMUTH, DEG. 

ELEVATION, DEG. 

LONGITUDE, DEG. 

LATITUDE, DEG. 

LONGITUE (RADAR) DEG 

0 a=-58. 15003437 

p=-. 00000113 

$=61 . 44889048 

0=18.44669380 

$5=77.00000000 

1 a=-57. 00645790 

P=3. 10616670 

$=64.07291615 

6=19.86970746 

$5=77.17186916 

2 a=-55. 34118427 

p=6. 99699965 

$=66.73900244 

6=21.25119328 

$s=77. 34373832 

3 a=-52. 71266671 

3=12.25898824 

<f>=69 . 44867889 

©=22.58706924 

*3=77.51560749 

4 a=-48. 00953651 

3=20.15528607 

$=72.20310392 

6=23.87327468 

$5=77.68747665 

5 a=-37. 56514360 

3=33.45800818 

$=75.00300642 

6=25.10579627 

$5=77.85934581 

6 a=-4. 22028134 

3=53.59066411 

$=77.84862649 

6=26.28069784 

$5=78.03121497 

7 a=64.33158348 

3=128.72017840 

$=80.73965749 

6=27.39415375 

$5=78.20308413 

8 a=90.55732702 

3=146.32232990 

$=83.67519231 

6=28.44248594 

$5=78.37495330 

9 a=99.36073597 

3=157.36021261 

$=86.65367667 

6=29.42220398 

$5=78.54682246 

10 0=103.52061532 

3=164.18076975 

$=89.67287325 

6=30.33004762 

$5=78.71869162 

11 0=105.90245379 

3=168.89052738 

$=92.72983999 

0=31.16303074 

$5=78.89056078 

12 0=107.42823026 

3=172.45952834 

$=95.82092612 

6=31.91848556 

$5=79.06242995 

13 0=108.47830717 

3=175.35471125 

$=98.94178883 

6=32.59410553 

$5=79.23429911 

14 0=109.23705338 

8=177.82162762 

$=102.08743269 

0=33.18798516 

$5=79.40616827 

15 0=109.80439856 

8=179.99999761 

$=105.25227275 

0=33.69865517 

$5=79.57803743 


Table 2. Representative data of entire package sent to iterator. 



(H, * 

o, i. Go, Ro ASSUMED KNOUN) 




T=122. 15382475 

0 o=-58. 15003437 

p=-. 00000113 



1=276.83607064 

1 o=-49. 49765192 

P=17. 81825854 



T=431. 51831654 

2 0=81.48848503 

P=41. 79462703 



1=586.20056243 

3 0=106.34412297 

p=169. 86594977 



1=740.88280832 

4 0=109.80439856 

P=179. 99999761 



Table 3. 

Results of iteration scheme. 


ORIGINRL GUESSES: 8=7000.00000 

kn Ecc=. 00000 Io=27.00000 

deg <*>=.00000 deg Q=. 00000 

deg 


ITERATION #0 R=6997.1192B kn Ecc=. 00153 <4=00.70312 deg Io=27. 63809 deg 0=359.66414 deg /[SUM (ERRORS) 2 ]=2009. 63630547 
ITERATION #1 R=6976. 37418 kn Ecc=. 01593 w=04. 53330 deg Io=34. 23148 deg 0=000.10470 deg /[SUM (ERRORS) 2 ]=1989. 44887953 
ITERATION #2 8=6984.14600 kn Ecc=. 02578 <4=04.30058 deg Io=34. 72687 deg 0=005.38340 deg /[SUM ( ERRORS) 2 ] =1784 . 30766470 
ITERATION 83 A=7001. 23236 kn Eec=. 03438 <4=05.26028 deg Io=35. 26693 deg 0=009.02111 deg /[SUM (ERRORS) 2 ]=1605. 87646748 
ITERATION 84 8=7023.14477 kn Ecc=. 04188 <4=06.72678 deg Io=35. 64320 deg 0=011.73783 deg /[SUM (ERR0RS) 2 ]=1445. 02945689 
ITERATION 85 8=7047.51661 kn Ecc=. 04841 <4=08.35303 deg Io=35. 86671 deg 0=013.91800 deg /[SUM (ERRORS) 2 ]=1300. 14138118 
ITERATION 86 8=7072. 8B632 kn Ecc=. 05411 <4=09.97658 deg Io=35. 97578 deg 0=015.75416 deg /[SUM (ERRORS) *]=1169. 70535455 
ITERATION 87 8=7098.31449 kn Ecc=. 05910 <4=11.51927 deg Io=36. 00527 deg 0=017.35192 deg /[SUM (ERRORS) 2 ]=1052. 31412419 
ITERRTION 88 8=7160.51885 kn Ecc=. 07010 <4=15.08418 deg Io=35. 94827 deg 0=020.90512 deg /[SUM (ERRORS) *]=8946. 67896590 
ITERATION 89 8=7212.08151 kn Ecc=. 07731 <4=17.84513 deg Io=35. 62406 deg 0=023.55637 deg /[SUM (ERRORS) 2 ] =0708 . 35808592 
ITERRTION 810 8=7254.11223 kn Ecc=. 08240 <4=19.81318 deg Io=35. 33142 deg 0=025.69427 deg /[SUM (ERRORS) 2 ]=8530. 02498390 
ITERATION 811 8=7288.44022 kn Ecc=. 08623 <4= 21.20240 deg Io=35. 13819 deg 0=027.41042 deg /[SUM (ERRORS) 2 ]=0396. 60902847 
ITERRTION 812 8=7316. 438BB kn Ecc=. 08921 <4=22.19628 deg Io=35. 02999 deg 0=028.75806 deg /[SUM (ERRORS) s ]=0296. 80609278 
ITERATION 813 8=7339.12253 kn Ecc=. 09158 <4=22.91851 deg Io=34. 97731 deg 0=029.79734 deg /[SUM (ERRORS) 2 ]=0222. 14316318 
ITERATION 814 8=7357.31780 kn Ecc=. 09347 <4=23.44940 deg 1 0=34.95647 deg 0=030.58963 deg /[SUM (ERRORS) 2 ]=0166.279115SS 
ITERRTION 815 8=7371. 751B7 kn Ecc=. 09497 <4=23.84247 deg I o=34 . 95223 deg 0=031.18950 deg /[SUM ( ERRORS) 2 ] =0124 .47395496 
ITERATION 816 8=7394.40775 kn Ecc=. 09732 <4=24.42709 deg Io=34. 95924 deg 0=032.09412 deg /[SUM (ERRORS) 2 ]=0093. 18528045 
ITERRTION 817 8=7406.78486 kn Ecc=. 09861 <4=24.71504 deg Io=34. 97549 deg 0=032.54899 deg /[SUM (ERRORS) 2 ]=0046. 30360731 
ITERRTION 818 R=7413. 29293 kn Ecc=. 09929 <4=24.85799 deg I 0=34.98695 deg 0=032.77596 deg /[SUM (ERRORS) 2 ]=0023. 02310822 
ITERRTION 819 8=7416.62822 kn Ecc=. 09964 <4=24.92917 deg Io=34. 99342 deg 0=032.88883 deg /[SUM ( ERRORS) 2 ]=0011 .45162358 
ITERRTION 820 8=7418.31322 kn Ecc=. 09982 <4=24.96465 deg Io=34. 99676 deg 0=032.94487 deg /[SUM (ERRORS) 2 ] =0005. 69708321 
ITERATION 821 8=7419.15825 kn Ecc=. 09991 <4=24.98236 deg I 0=34. 99843 deg 0=032.97267 deg /[SUM (ERRORS) 2 ] =0002 . 83453385 
ITERATION 822 A=7419.5B049 kn Ecc=. 89995 <4=24.99119 deg Io=34. 99924 deg 0=032.98645 deg /[SUM (ERRORS) 2 ]=0001 .41037319 
ITERRTION 823 8=7420.00172 kn Ecc=. 10000 <4=25.00001 deg Io=35. 00003 deg 0=033.00012 deg /[SUM (ERRORS) 2 ]=0000 .70177748 
ITERATION 824 8=7420.00025 kn Ecc=. 10000 <4=25.00000 deg Io=35. 00000 deg 0=033.00000 deg /[SUM (ERRORS) 2 ] =0000.00353102 
ITERRTION 825 8=7420.00026 kn Ecc=. 10000 <4=25.00000 deg Io=35. 00000 deg 0=033.00000 deg /[SUM ( ERRORS) 2 ] =0000 . 80001940 
ITERRTION 826 8=7420.00026 kn Ecc=. 10000 <4=25.00000 deg Io=35. 00800 deg 0=033.00000 deg /[SUM (ERRORS) 2 ]=0000. 00000010 


CONVERGED SOLUTION: 8=7420.00026 kn Ecc=. 10000 Io=35. 00000 deg <4=25.00000 deg 0=33.00000 deg 

/[SUM (ERRORS) 2 ]=. 00008010 
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Table 4. Second example. 


■MAJOR AXIS (A)=6019. 

KM, ECC=. 72428 INC(I)=.000 

DEG., ARG. PERIGEE(u) 

=225.000 DEG., LONG. ASC. 

MODE (fi) = . 000 DEG. 

APPARITION OCCURS AT 

f =154 . 991 deg. [T=1128.749 

sec.] AND OCCULTATION 

OCCURS AT f =219 . 128 deg. 

[T =3932. 605 sec.] 

TRUE ANOMOLY, DEG. 

RADIUS, KM. 

RANGE, KM. 

ECC. ANOMOLY, DEG. 

TIME, SEC. 

0 f =154. 99056464 

r=8327. 03917718 

P=5353. 56799812 

£=121 .97277140 

T=1120. 74879014 

1 f=159. 84812332 

r=8940. 30401130 

P=5316. 690801 46 

£=132.08228159 

1=1308.20507948 

2 f=164. 13710554 

r=9434. 24059567 

p=5301. 21356239 

£=141.58258033 

T=1495. 66296882 

3 f=168. 04340171 

r=9810.4403B207 

p=5278. 16143365 

£=150.64963478 

T=1683. 12005817 

4 f =171. 69290062 

r=10099. 63885244 

p=5226. 61914198 

£=159.41388293 

T=1870. 57714751 

5 f=175. 17806710 

r=l 0282. 4 83021 65 

P=5131. 41085411 

£=167.97870649 

1=2058.03423685 

6 f=178. 57269911 

r=10369. 54893382 

P=4981. 20243176 

£=176.43164170 

T=2245. 49132620 

7 f=181. 94118337 

r=10362. 37915844 

P=4767. 10213407 

£=184.85198369 

T=2432. 94841554 

8 f=185. 34530128 

r=10260. 78748541 

p=4481. 69994810 

£=193.31671190 

T=2620. 40550488 

9 f=188. 85039710 

r=18063. 25051039 

P=4118. 50129078 

£=201.90596867 

T =2807. 86259422 

10 f =192. 53242283 

r=9766. 70455542 

p=3671. 85978114 

£=210.70915095 

T=2995. 31968357 

11 f=196. 48778574 

r=93b6. 28463255 

P=3137. 99737918 

£=219.83291254 

T =3182.77677291 

12 f=200. 84937183 

r=8B54. 87299320 

P=2519. 60614320 

£=229.41318022 

1=3370.23386225 

13 f=205. 81589808 

r=8222. 35862730 

P=1B46. 25538555 

£=239.63526416 

T=3557. 69095160 

14 f=211. 71283239 

r=7454. 43313420 

p=1281. 26945377 

£=250.77108362 

1=3745.14804094 

15 f= 219.12843774 

r=6530. 64834910 

p=1403. 73789664 

£=263.25619535 

T=3932. 60513028 

AZIMUTH, DEG. 

ELEURTION, DEG. 

LONGITUDE, DEG. 

LATITUDE, DEG. 

L0NGI1UE (RADAR) DEG 

0 c=-90. 00000000 

3=. 00000600 

4>=19 . 99056464 

©=.00000000 

<{>3=68.00000000 

1 a=“90. 00000000 

3=9.31971328 

<{>=24.84812332 

©=.00000000 

$3=60.78107121 

2 a=-90. 00000000 

3=17.40117751 

4=29.13710554 

©=.00000000 

$s=61. 56214241 

3 a=-90. 00000000 

3=24.44702159 

4=33.04340171 

©=.00000000 

$5=62.34321362 

4 a=-90. 00000000 

3=30.66801485 

4=36.69290062 

©=.00000000 

$3=63.12428482 

5 a=-90. 00000000 

3=36.26370188 

4=40.17806710 

©=.00000000 

$s=63.90535603 

6 a=-90. 00000000 

3=41.41969757 

4=43.57269911 

©=.00000000 

$3=64.68642723 

7 a=-90. 00000000 

3=46.31600507 

4=46.94118337 

©=.00000000 

$3=65.46749844 

8 a=-90. 00000000 

3=51.14502534 

4=50.34530128 

©=.00000000 

$3=66.24856964 

9 a=-90. 00000000 

8=56.14476239 

4=53.85039710 

©=.00000000 

$3=67.02964085 

10 a=-90. 00000000 

3=61.66664772 

4=57.53242283 

©=.00000000 

$3=67.81071206 

11 a=-90. 00000000 

3=68.33790782 

4=61.48778574 

©=.00000000 

$3=68.59178326 

12 o=-90. 00000000 

3=77.52662540 

4=65.84937183 

©=.00000000 

$s=69 . 37285447 

13 a=-90. 00000000 

3=87.05064752 

4=70.81589808 

©=.00000000 

$3=70.15392567 

14 a=-90. 00000000 

3=125.84723511 

4=76.71203239 

©=.00000000 

$3=70.93499688 

15 a=-90. 00000000 

3=179.99998791 

4=84.12843774 

©=.00000000 

$3=71.71606808 


Table 5. Representative data of entire package sent to iterator. 

(H, *, 00, Ro ASSUMED KNOUM) 

T=1120 . 74B79014 p=. 00000600 

1=2526.67696021 0=48.72566711 

1=3932.60513028 0=179.99998791 

Table 6. Results of iteration scheme. 

ORIGIMRL GUESSES: A=637B. 00008 kn Ecc=. 00000 Io=. 00000 deg u=. 00000 deg 0= . 00000 deg 


ITERATION 

80 

ITERATION 

81 

ITERATION 

«2 

ITERATION 

#3 

ITERATION 

84 

ITERATION 

#5 

ITERATION 

#6 

ITERATION 

87 

ITERATION 

88 

ITERATION 

#9 

ITERATION 

#10 

ITERATION 

#11 

ITERATION 

#12 

ITERATION 

#13 

ITERATION 

#14 

ITERATION 

#15 

ITERATION 

#16 

ITERATION 

#17 

ITERATION 

#18 

ITERATION 

#19 

ITERATION 

#20 

ITERATION 

#21 

ITERATION 

#22 

ITERATION 

#23 

ITERATION 

#24 

ITERATION 

#25 

ITERATION 

#26 

ITERATION 

#27 

ITERATION 

#28 

ITERATION 

#29 

ITERATION 

#30 

ITERATION 

831 

ITERATION 

#32 

ITERATION 

#33 

ITERATION 

#34 

ITERATION 

#35 


0=6377.82547 kn Ecc= 
A=6377. 65093 kn Ecc= 
R=6377. 47640 kn Ecc= 
R=6377. 30187 kn Ecc= 
R=6377. 47640 kn Ecc= 
R=6377. 30187 kn Ecc= 
R=6377. 47640 kn Ecc= 
R=6377. 30187 kn Ecc= 
R=6377. 47640 kn Ecc= 
R=6377. 65093 kn Ecc= 
R=6377. 82547 kn Ecc= 
fl=6377. 65093 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6377. 65093 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=637B. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=6377. 82547 kn Ecc= 
R=6378. 00000 kn Ecc= 
R=637B. 00000 kn Ecc= 
R=6378. 00000 kn Ecc= 


.01000 w=246. 53286 
.11000 u=275. 73814 
.21000 w=243. 58550 
.31000 U=214. 12491 
.41000 14=258.37404 
.51000 (4=119.31538 
.41000 (4=161.91792 
.31000 (4=143.56112 
.21000 (4=002.73382 
.11000 (4=102.47098 
.01000 (4=218.38686 
.11000 (4=076.02424 
.01000 (4=060.90440 
.11000 (4=110.29098 
.01000 (4=133.66279 
.01000 (4=110.56111 
.01000 (4=220.70121 
.11000 (4=336.81249 
.01000 (4=056.21260 
.11000 (4=036.52690 
.01000 (4=107.39903 
.01000 (4=238.21842 
.11000 <4=097.93358 
.01000 (4=262.73277 
.01000 (4=037.54403 
.01000 <4=074.54306 
.01000 (4=274.77163 
.01000 <4=0 47.80174 
.01000 (4=322.64812 
.01000 (4=062.38489 
.11000 (4=116.39484 
.01000 (4=059.07265 
.11000 (4=354.23826 
.01000 (4=326.15469 
.01000 <4=028.42436 
.01000 (4=252.93308 


deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 
deg /[SUM (ERRORS) 


*]=3342. 25372390 

•]=1707. 90104467 

•]=1563. 56658065 

•J =0967 .94274041 

•]=1218. 41028906 

8 ] =0862 .33919423 

2 ]=4237. 85560753 

•] =3587. 491 671 98 

•]=3902. 1231 3582 

8 ]=3642. 75359697 

8 ]=3958. 24350289 

•3=1850.50353792 

•]=4175. 96923678 

•3=5569.74171210 

•3=3903.96875415 

•3=3621.66578190 

•3=3759.25058461 

•3=1828.42076689 

•3=3104.59657251 

•3=4543.88574404 

•3=4167.53932133 

•3=3783.94734648 

•3=1756.34696722 

•3=3991.99014902 

•3=1758.09479345 

•3=4138.43426753 

•3=4115.06597338 

•3=1855.01818369 

•3=4360.17724804 

•3=2705.93322260 

•3=5573.56804375 

•3=3864.88321210 

•3=4606.71200290 

•3=3371.78846483 

•3=2762.87190362 

•3=3943.68348011 
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Table 6. Results of iteration scheme (continued). 


ITERATION 836: R=6378. 08000 kn Ecc=. 01000 <4=064.59755 deg /[SUN (ERRORS) *3=1715.94468123 

ITERRTION 837: 8=6377.82547 kn Ecc=. 11000 <4=125.01488 deg /[SUM (ERRORS) *3=5579. 52089429 

ITERATION 838: 8=6378.00000 kn Ecc=. 01000 (4=295.14937 deg /[SUM (ERRORS) *3=3814.22587824 

ITERATION 839: R=6378. 00000 kn Ecc=. 01000 (4=059.07106 deg /[SUM (ERRORS) *3=2112.08775784 

ITERATION 840: R=6377. 82547 kn Ecc=. 11000 <4=354.26307 deg /[SUM (ERRORS) *3=4686.67699080 

ITERATION 841: 8=6378.00000 kn Ecc=. 01000 <4=325.54414 deg /[SUM (ERRORS) *3=3372.17423267 

ITERRTION 842: 8=6378.00000 kn Ecc=. 01000 <4=034.50583 deg /[SUM (ERR0RS)*]=2752. 92595544 

ITERATION 843: R=6378. 00000 kn Ecc=. 01000 <4=132.14604 deg /[SUM (ERRORS) *3=4073 .28820556 

ITERRTION 844: R=6378. 00000 kn Ecc=. 01000 <4=121.06298 deg /[SUM (ERRORS) *3 =3628. 40253180 

ITERATION 845: R=6378. 00000 kn Ecc=. 01000 <4=181.18194 deg /[SUM (ERRORS) *3=3687 .17028436 

ITERRTION 846: R=6377. 82547 kn Ecc=. 11000 <4=094.81280 deg /[SUM (ERRORS) *3=3172.60797744 

ITERRTION 847: 8=6378.00000 kn Ecc=. 01000 <4=291.13616 deg /[SUM (ERRORS)*3=4015. 96918752 

ITERRTION 848: R=6378. 00000 kn Ecc=. 01000 <4=146.17872 deg /[SUM (ERRORS) *3=2053.68258407 

ITERATION 849: 0=6378.00000 kn Ecc=. 01000 <4=343.47329 deg /[SUM (ERRORS)*3=3580. 81284395 

ITERRTION 858: R=6378. 00000 kn Ecc=. 01000 <4=218.86757 deg /[SUM (ERRORS) *3=3030 . 67410247 

ITERRTION 851: 8=6377.82547 kn Ecc=. 11000 <4=055.29646 deg /[SUM (ERRORS) *3=1845 . 65847847 

ITERRTION 852: 0=6378.00000 kn Ecc=. 21000 <4=048.84146 deg /[SUM (ERRORS) *3=4523 .50960091 

ITERRTION 853: 8=6377.82547 kn Ecc=. 31000 <4=260.19530 deg /[SUM (ERRORS) *3=4415.62354477 

ITERATION 854: R=6377. 65093 kn Ecc=. 41000 <4=038.99094 deg /[SUM (ERRORS) *3=0964 . 67445113 

ITERRTION 855: R=63?7. 47640 kn Ecc=. 31000 <4=211.81071 deg /[SUM (ERRORS) *3=4293.98179749 

ITERRTION 856: R=6377. 65093 kn Ecc=. 41000 <4=263.74250 deg /[SUM (ERRORS) *3=1283.33023873 

ITERRTION 857: 8=6377.47640 kn Ecc=. 51000 <4=083.08813 deg /[SUM (ERRORS) *3=1011 .94486805 

ITERATION 858: 8=6377.30187 kn Ecc=. 61000 <4=312.47459 deg /[SUM (ERRORS) *3 =4456 .22739218 

ITERRTION 859: 8=6377.12734 kn Ecc=. 71000 <4=138.56418 deg /[SUM (ERRORS) *3=2377 .60300719 

ITERRTION 860: R=6377. 30187 kn Ecc=. 81000 <4=213.83595 deg /[SUM (ERRORS) *3 =3887 .18052500 
ITERRTION 861: R=6377. 47640 kn Ecc=. 75529 <4=225.56000 deg /[SUM (ERRORS)*3=1250. 86099995 
ITERATION 862: 8=6277.47640 kn Ecc=. 71739 <4=224.33091 deg /[SUM (ERRORS) *3=8460 . 12395739 
ITERRTION 863: 8=6177.47640 kn Ecc=. 71899 <4=224.13317 deg /[SUM (ERRORS) *3 =0330.28472467 
ITERATION 864: 8=6077.47640 kn Ecc=. 72094 <4=224.66098 deg /[SUM (ERRORS) *3=0217.45002740 
ITERATION 865: 8=6002.18338 kn Ecc=. 72299 <4=224.86027 deg /[SUM (ERRORS) *3 =8080 . 15408107 
ITERATION 866: 8=6022.73748 kn Ecc=. 72463 <4=225.13766 deg /[SUM ( ERRORS) *3 =0020 . 36231904 
ITERRTION 867: R=6016. 76933 kn Ecc=. 72419 <4=224.90929 deg /[SUM (ERRORS) *3=0000.57612011 
ITERRTION 868: 8=6020.10635 kn Ecc=. 72433 <4=225.05943 deg /[SUM ( ERRORS) *3 =0085 .34985398 
ITERRTION 869: 8=6017.80717 kn Ecc=. 72425 <4=224.96052 deg /[SUM (ERRORS) *3=0083. 48845386 
ITERRTION 870: 8=6019.36574 kn Ecc=. 72430 <4=225.02634 deg /[SUM ( ERRORS) *3 =0002 .32481 090 
ITERRTION 871: R=601B. 32287 kn Ecc=. 72427 <4=224.98242 deg /[SUN (ERRORS) *3=0001 .55106925 
ITERRTION 872: 8=6019.01815 kn Ecc=. 72429 <4=225.01173 deg /[SUM (ERRORS) *3=0001 .03537443 
ITERATION 873: R=601B. 55427 kn Ecc=. 72427 <4=224.99218 deg /[SUM (ERRORS) *3 =0000.69857852 
ITERRTION 874: 8=6018.86368 kn Ecc=. 72429 <4=225.00522 deg /[SUM (ERRORS) *3=0000. 46077383 
ITERATION 875: 8=6018.65724 kn Ecc=. 72428 <4=224.99652 deg /[SUM (ERRORS) *3=0000.30736057 
ITERATION 876: 8=6018. 79495 kn Ecc=. 72428 <4=225.00232 deg /[SUM (ERRORS) *3=0000.20586258 
ITERRTION 877: 8=6018.70308 kn Ecc=. 72428 <4=224.99845 deg /[SUM (ERRORS) *3=0000. 13679601 
ITERRTION 878: 8=6018.76437 kn Ecc=. 72428 <4=225.00103 deg /[SUM (ERRORS) *3=0008.09126299 
ITERRTION 879: 8=6018.72348 kn Ecc=. 72428 <4=224.99931 deg /[SUM (ERRORS) *3=0000.06888258 
ITERRTION 880: 8=6018.75075 kn Ecc=. 72428 <4=225.00046 deg /[SUM (ERRORS) *3=0000.04061689 
ITERRTION 881: 8=6018.73256 kn Ecc=. 72428 <4=224.99969 deg /[SUM (ERRORS) *3=0000.02709631 
ITERRTION 882: 8=6018. 74470 kn Ecc=. 72428 <4=225.00021 deg /[SUM (ERRORS) *3=0000 .01887675 
ITERRTION 803: R=6018. 73660 kn Eec=. 72428 <4=224.99986 deg /[SUM (ERRORS) *3=0000.01205941 
ITERATION 884: R=6018. 74200 kn Ecc=. 72428 <4=225.00089 deg /[SUM (ERRORS) *3=0000.00804516 
ITERRTION 885: R=6018. 73840 kn Ecc=. 72428 <4=224.99994 deg /[SUM (ERRORS) *3=0000.08536712 
ITERATION 806: R=6018. 74080 kn Ecc=. 72428 <4=225.00004 deg /[SUM (ERRORS) *3=0000.00358055 
ITERRTION 887: R=6018. 73920 kn Ecc=. 72428 0=224.99997 deg /[SUM (ERRORS) *3=0000.00238867 
ITERRTION 888: R=6818. 74027 kn Ecc=. 72428 <4=225.00002 deg /[SUM (ERRORS) *3=0000.00159355 
ITERRTION 889: R=6018. 73955 kn Ecc=. 72428 <4=224.99999 deg /[SUM (ERRORS) *3=0000.00106309 
ITERRTION 890: 8=6018.74003 kn Ecc=. 72428 <4=225.00001 deg /[SUM (ERRORS) *3=0000.00070922 
ITERRTION 891: R=6018. 73971 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3 =0000 .00847314 
ITERRTION 892: R=6018. 73992 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00031564 
ITERRTION 893: R=6018. 73978 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00021057 
ITERRTION 894: 8=6018.73988 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3 =0000 .00814048 
ITERRTION 895: R=6010. 73981 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00009372 
ITERRTION 896: 8=6018.73986 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00006252 
ITERATION 897: 8=6018.73983 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3 =0000 .00084171 
ITERRTION 898: 8=6018.73985 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00002783 
ITERRTION 899: 8=6018.73983 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000 .00001856 
ITERRTION 8100: R=601B. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00001238 
ITERRTION 8101: 8=6018.73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00000826 
ITERRTION 8102: 8=6018.73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS)*3=0000. 00000551 
ITERRTION 8103: R=6018. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS)*3=0000.0000036'fe 
ITERRTION 8104: R=6018. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS)*3=0B00. 00000245 
ITERRTION 8105: R=6018. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3 =0000 .00000164 
ITERRTION 8106: R=6018. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3=0000.00000109 
ITERATION 8107: R=6018. 73984 kn Ecc=. 72428 <4=225.00000 deg /[SUM (ERRORS) *3 =0000. 00000073 


CONVERGED SOLUTION: fl=6018. 73984 kn Ecc=. 72428 I o=. 00000 deg <4=225.00000 deg fi=. 00000 deg 

/[SUM (ERRORS) *3=. 00000073 
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Figure 1. Coordinate systems showing azimuth (a) and elevation (/?). 
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Figure 2. Orbital elements. 
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